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Abstract 

We complete the non-perturbative calculations of the strange quark mass 
and the Lambda parameter in two flavor QCD by the ALPHA collaboration. 
The missing lattice scale is determined via the kaon decay constant, for whose 
chiral extrapolation complementary strategies are compared. We also give a value 
for the scale r in physical units as well as an improved determination of the 
renormalization constant Z\. 
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1 Introduction 



The parameters of the standard model of particle physics have to be determined 
by matching theory to experimental data. For two flavor QCD we present results 
for two of these parameters, the scale parameter Aq CD and the mass of the strange 
quark. 

These results are the outcome of a long project by the ALPHA collabora- 
tion. The general strategy using the Schrodinger functional to define the coupling 
constant has been laid out in Refs. [1,2], and results for the A parameter in pure 
gauge theory have been published in [3]. In the two flavor theory for A [4,5] and 
the strange quark mass [6] , however, the determination of the lattice scale, which 
allows for their conversion to physical units, has been lacking. This determination 
is the subject of the present paper. 

In lattice computations, the physical mass scale is set by picking one dimen- 
sionful observable and identifying its value at the point where the quark mass 
ratios correspond to the physical situation with the experimental input. If we 
have done the calculation with all physical effects taken into account (and if the 
theory is correct), it does not matter, which observable we take. Here, however, 
we restrict ourselves to QCD with two dynamical flavors of light quarks leading 
to a systematic uncertainty which is hard to determine. In our computation, we 
use the kaon decay constant to set the scale. Over the pion decay constant it 
has the advantage of a chiral extrapolation which is milder and therefore better 
under control. 1 However, we need a quenched strange quark. Also the mass of 
the Omega baryon is popular to set the scale and first results indicate that similar 
numbers are obtained from this observable [7]. 

In previous publications the results have been converted to physical units 
using the scale parameter r [8], defined via the force between static quarks. The 
conversion relied on measurements of r /a by QCDSF [9] and the assumption that 
r = 0.5 fm. Below, we present our own results for r /a, which differ substantially 
from the previous values and lead to an update in r A. 

The paper is organized as follows. The lattice action, an overview of the 
ensembles and details of the error analysis are given in Sect. 2, followed by the 
definition of the hadronic observables in Sect. 3 and results for the scale parameter 
r in Sect. 4. The strategies for the chiral extrapolation are discussed in Sect. 5 
leading to the scale determination from the kaon decay constant. The results for 
the Lambda parameter and the strange quark mass are contained in Sects. 6 and 
7, respectively. 

1 This property depends on how one actually approaches the physical point. It is in particular 
true for the strategy 1 which we introduce below. 
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The appendices contain updates of many quantities, whose analysis has been 
subject of previous publications. The renormalization constants and Zp are 
discussed in App. B and App. C, the hadronic scale L\ of the Schrodinger func- 
tional calculations is subject of App. D, followed by a determination of the critical 
mass of the improved Wilson fermions and the singlet renormalization factor in 
App. E. 

2 Lattice parameters and simulation algorithms 

In this computation we use the Wilson plaquette gauge action for the gluon fields 
together with two degenerate flavors of 0(a) improved Wilson fermions [10]. The 
action 

S[U,4,il>] = (3^ti[l-U{p)}+a A ^P{x){D + m Q )4>{x) (2.1) 

p X 

has three parameters: g , m and c sw . The coupling constant g is given by 
(3 = 6/<7q. The fermions with bare mass m , usually substituted by the hopping 
parameter k — (8 + 2am )~ 1 , are implemented by the lattice Dirac operator 

D = ^(V; + V„) - aV;V,} + c SWJ a, u F, u . (2.2) 

It includes the covariant forward and backward derivatives, and V* , and the 
Sheikholeslami-Wohlert [1 1] improvement term involving the standard discretiza- 
tion F^ u of the field strength tensor [12]. Its coefficient c sw has been determined 
non-perturbatively [13]. 

We have generated ensembles at three values of f3 — 5.2, 5.3 and 5.5 which 
correspond roughly to lattice spacings of a = 0.076 fm, 0.066 fm and 0.049 fm, 
respectively, with details given in Section 5.6. The ensembles are listed in Table 1. 
All lattices have size (2L) x L 3 and the pion mass is always large enough such 
that ffljjL > 4. We therefore expect finite size effects to be small. 

2.1 Simulation algorithms 

For most of the ensembles, generated within the CLS effort, 2 the DD-HMC 
algorithm [14, 15] has been used as implemented in the software package by 
M. Liischer [16]. It is based on a domain decomposition to separate the infrared 
from the ultraviolet modes of the fermion determinant. A main feature is the 
locally deflated, Schwarz preconditioned GCR solver [17, 18] which significantly 
reduces the increase in computational cost as the quark mass is lowered. 

2 https : //twiki . cern . ch/twiki/bin/view/CLS/ 
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id 


L/a 


/3 


K 


K s 


Rq 


m^MeV] 




A2 


32 


5 2 


1 3565 


1 35438(20"! 


5 485(21 1 


630 


7 7 


A3 






0.13580 


0.135346(20) 


5.674(32) 


490 


6.0 


A4 






0.13590 


0.135285(20) 


5.808(34) 


380 


4.7 


A5 






0.13594 


0.135257(20) 


5.900(24) 


330 


4.0 


E4 


32 


5.3 


13610 


1 35836 ( 17) 




580 


6.2 


E5 






0.13625 


0.135777(17) 


6.747(59) 


440 


4.7 


F6 


48 




0.13635 


0.135741(17) 


6.984(51) 


310 


5.0 


F7 






0.13638 


0.135730(17) 


7.051(43) 


270 


4.3 


N4 


48 


5.5 


0.13650 


0.136278(08) 


9.32(30) 


550 


6.5 


N5 






0.13660 


0.136262(08) 


9.31(26) 


440 


5.2 


N6 






0.13667 


0.136250(08) 


9.55(11) 


340 


4.0 


07 


64 




0.13671 


0.136243(08) 


9.68(10) 


270 


4.2 



Table V. Overview of the ensembles used in this study. We give the label, the spatial extent 
of the lattice, /3 — 6/g$, the hopping parameter k of the sea quarks, the hopping parameter k s 
of the strange quark, the scale i? = r /a, the mass of the sea pion and the product to^L, 
which is always larger than or equal to 4. All lattices have dimension T x L 3 with T = 2L. 

The drawback of this algorithm is that due to the block decomposition only a 
fraction of gauge links -Ractive is updated during a trajectory. In pure gauge theory 
the autocorrelation times are inversely proportional to this fraction of active links 
[19]; we expect this behavior also in the theory with fermions. Typical domain 
decompositions lead to active link ratios between 0.37 and 0.5 and therefore a 
factor between 2 and 3 increased autocorrelation times. 

For some lattices, we therefore used a Hybrid Monte Carlo algorithm [20] with 
a mass preconditioned fermion determinant [21,22]. Our implementation [23], MP- 
HMC from here on, is based on the DD-HMC package and in particular takes over 
the deflated solver because of its efficient light quark inversions. This algorithm 
was employed for ensembles A5, N6 and 07 given in Table 1; all other ensem- 
bles were generated with the DD-HMC. Appendix A gives details about these 
algorithms and the values of the parameters used for the gauge field generation. 

2.2 Autocorrelations 

In Monte Carlo data the effect of the autocorrelations has to be accounted for in 
the error analysis. For all observables F = F(a 1: . . . , a n ), functions of expectation 
values dj = (A) of primary observables Aj, we therefore compute an estimator of 
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the autocorrelation function 



rV(0 = E /i(a)/i(3)<(A(*) - Oi)(^-(0) - a,)) (2.3) 



where /j = diF(a) following the procedures detailed in Ref. [24]. The argument t 
indicates the Monte Carlo time. The integrated autocorrelation time is then 



(F) = - + Y^t\- (2-4) 



Tint 

Z 

t=l 



which then enters the statistical error of the observable op from N measurements 

4 = 2^IV(0) . (2.5) 

The sum in eq. (2.4) is normally truncated at a "window" W [25] which bal- 
ances the statistical uncertainty due to the limited sample size and the systematic 
error coming from neglecting the tail for t > W. The value of W is determined 
from the measurement of pF{t) = IV(t)/rV(0) alone and for each F separately. 
Neglecting the tail above W leads — at least on average — to an underestimation 
of Tint and the statistical error of the observable. It is particularly problematic 
in the presence of slow modes of the Monte Carlo transition matrix which only 
couple weakly to the observable in question. To account for them we use the 
method outlined in Ref. [19], estimating their time constants from observables to 
which the slow modes couple strongly. Using them, we can then estimate the tails 
of the autocorrelation functions of the observables we are interested in and arrive 
at a more conservative error estimate. 

Experience tells us that for small lattice spacing the topological charge is 
particularly sensitive to slow modes [19,26], for which we use the field theoretical 
definition after smoothing the field by the Wilson flow integrated up to to as 
defined in Ref. [27]. Actually, only the square of the charge needs to be considered, 
because we are only interested in parity even observables. Unfortunately we are 
not in the position to accurately determine its autocorrelation time for most of 
our ensembles. We therefore combine the scaling laws found in pure gauge theory 
[19] with the measurement for our high statistics ensembles E5 and arrive at the 
estimate 

e 7(/3-5.5) 

r cxp (/3) = 200^- , (2.6) 

-'•'active 

in units of molecular dynamics time with c T = 2 for trajectories of length r = 0.5 
and c T = 1 for r = 2 and 4. The values of inactive can be found in Table 9 for the 
DD-HMC algorithm and is equal to one for the MP-HMC. 
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An example of the procedure is given in Fig. 1, showing the autocorrelation 
function of the kaon decay constant Fk on the 07 ensemble. Using the standard 
procedure [24,25], the sum in eq. (2.4) is truncated at the window Wi from which 
we would get r; nt = 0.7. When the contribution of the tail is included, the im- 
proved estimate gives r int = 4, which translates to a more than doubled error 
estimate. 




W u Wy 



200 



Figure 1: Autocorrelation function of Fk for the 07 lattice. The line gives our estimate for 
its tail. The standard method of Ref. [24] gives a window W = Wi and r; n t = 0.7, compared to 
Tint = 4 including the tail contribution which we add from W = W u , or more than a factor two 
in the error of the observable. 



3 Observables 

The kaon decay constant necessarily requires the addition of a quenched strange 
quark to the Nf = 2 theory. We denote the hopping parameter of this third flavor 
as k 3 and have k% = n 2 = K sea for the two sea quarks. For the purpose of a 
definition of the strange quark mass from the PCAC relation, we do in fact add a 
fourth quenched flavor with k 4 = k 3 . 

The computation of our pseudoscalar observables is based on two-point func- 
tions of the pseudoscalar density and the time component of the axialvector cur- 
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rent. At a fixed /t sca they are constructed from two valence quarks r and s 

m^o) = -« 3 X> rS (^ Sr (0)> (3-1) 

X 

/X^o) = -a 3 X>oW sr (0)> (3.2) 

X 

with P rs = ipr'-fsips and A™ = VVToTsVv This notation and the analysis that 
follows is similar to the one presented in Ref. [28]. Using the PCAC relation, 
average quark masses of flavors r and s can then be defined as 3 

, x U d o + d o)fAp(x ) + c A ad^d fp P (x ) 
m rs (x ) = ^-^ . (3.3) 

In this formula, do and <9q denote the forward and backward difference operators 
in time direction. The improvement coefficient c A has been determined non- 
perturbatively [29]. 

For sufficiently large x the mass m rs (x Q ) will have a plateau over which we 
can average. From its value m rs the renormalized quark mass is obtained [30] 

rs _ Z A (l + 6 A am se a + b A am rs ) 
m K - - — ~ -m rs , (d.4J 

Z P (1 + bpam sca + bpam rs ) 

with m sca = mi2- We will use one- loop perturbation theory for the improvement 
coefficients b A , bp, b A and bp, noting that they multiply very small terms. At 
this order in perturbation theory b A = bp = 0, while b A = 1 + 0.06167 g$ and 
bp = 1 + 0.06261 gl computed from the perturbative coefficients of [31]. An 
update of the non-perturbative determination of Z A [32] and Zp [6] is given in 
App. B and App. C, respectively. 

The renormalization and improvement of the PCAC quark masses is much 
simpler than the corresponding expression in terms of the bare subtracted quark 
masses m^ r = mo >r — m cr , where terms proportional to m sca are present already at 
the leading order in a [30]. In our analysis we therefore only use the renormalized 
PCAC relation eq. (3.4). The alternative definition of renormalized quark masses 
as well as the determination of additive renormalization m CT and the multiplicative 
renormalization factor Z m is discussed in App. E. 

3.1 Computation of the two-point functions 

We compute the two-point functions eq. (3.2) using U(l) noise sources f]t{x) = 
5 t>xo exp(i<f>(x)) located on randomly chosen time slices t [33,34]. Solving the Dirac 



3 This definition differs by a factor of two from Ref. [28] . 
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equation once for each noise vector Q = Q~ 1 (m Qr )ri t = a~ l (D + mo, r ) _1 75?7t 
is sufficient to get an estimator for the two-point functions projected to zero 
momentum 

a 3 /ppK) = J2^t(^ + t,xfC t s (x + t,x)), (3.5) 

X 

« 3 /Xp(^o) = J2^t^o + t,S)HQ(x Q + t,x)) (3.6) 

x 

where the average is over noise sources and gauge configurations. For our lattices, 
we use 10 noise sources per configuration, balancing the numerical cost and the 
accuracy which we wanted to reach on the given ensembles. 

3.2 Analysis of the data 

The following presentation applies to any flavor combination "rs" and we drop 
this sub/superscript for the sake of brevity. The mass of the pseudoscalar meson 
mps and its decay constant /ps can be extracted from /pp and the PCAC mass. 
For infinite time extent T, the spectral decomposition gives an expansion in terms 
of functions which decay exponentially for large time separations 

oo 

/pp(xo) = J>e-^° (3.7) 
i=i 

with E 1 = m PS the energy of the ground state and E 2 < E 3 < . . . the excited 
state contribution. For large time separations, we can then extract the decay 
constant from the leading coefficient 

f PS = Z A (1 + b A am sca + ab A m rs ) fps c , (3.8) 
f^ e = 2v^m rs m re 3/2 . (3.9) 

In the analysis with finite time extent T and time separation x we have 
to deal with particles running backwards in time and excited states. Since our 
lattices are large and the statistical precision of pseudoscalar correlators does not 
deteriorate at large x , we use the following procedure to fix the region x Q G 
[x™ in ,T — £™ in ], in which we can neglect the excited state contribution: we first 
perform a fit to the data using the first two terms in the expansion eq. (3.7), now 
including the finite T effects 

/pp(xo) = ci [e- ElX » + e ~ E ^ T ^] + c 2 [e~ E ^ + e -^( T -*o)] ( 3 .1 ) 

to a range where this function describes the data well, given the accuracy of the 
data. We then determine x™ m : it is the smallest value of x where the statistical 
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Xo/a 

Figure 2: The effective pion mass given by cosh(M cff (i - T/2))/ cosh(Af cff (t + 1 - T/2)) = 
fpp(t)/fpp(t + 1) for the 07 lattice. A two-state fit to data outside the shaded area determines 
^min rp ne regu j|. £ ^ ne g na i one-state fit is given by the error band. 

uncertainty on the effective mass is four times larger than the contribution of 
the excited state as given by the result of the fit using eq. (3.10). In a second 
step, only the first term of eq. (3.10) is fitted to the data restricted to this region. 
Figure 2 illustrates the procedure on our largest lattice. Formally at large xq 
and small sea quark mass the leading correction to eq. (3.7) comes from states 
which additionally to the ground state have two pions, E 2 ~ E\ + 2m n . However, 
for small quark masses and large L, the coefficient c 2 can be computed in chiral 
perturbation theory. It turns out to be very small [35] (at least for our large 
volumes) and such a contribution is invisible within our precision. In this sense 
the value of E 2 determined by the fit may actually be a higher state, which is one 
reason for us to use this excited state fit only in order to determine a safe x™ m . 

The error on the decay constant and the meson mass mps = E\ are computed 
including the autocorrelations as described in Sec. 2.2. In appendix F we give the 
partially quenched meson masses, decay constants and quark masses as well as 
the values of x™ in , which we use in our analysis. 

4 Scale parameter rg 

The analysis strategy for the scale Rq(/3,k) = r^/a [8] is based on [36] and we 
refer to this work for more detailed explanations and notation. The procedure con- 
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(3 n\ n 2 n 3 



5.2 
5.3 
5.5 



6 
8 
16 



9 

12 
24 



15 
20 
40 



Table 2: The number of smearing levels n\ used to construct the Operators Oi, I = 1,2,3 for 
measuring the Wilson loop correlation matrix. 

sists of measuring on-axis Wilson loops W(r,t) on smeared gauge configurations, 
extracting the static potential V(r) and finally solving the equation 



where F(r) = V'(r) is the static force. For the latter we use an improved definition 
which eliminates cut-off effects at tree level. Wilson loops are measured on gauge 
configurations after all links are replaced by HYP smeared [37] links. We take the 
HYP2 parameter choice [38] a x = 1.0, a 2 = 1.0 and a 3 = 0.5 and do one HYP- 
smearing level. The Wilson loops can be exactly represented as an observable in 
a theory including static quarks and this first smearing step corresponds to the 
choice of the static quark action (as far as the time-like links are concerned). On 
the HYP2-smeared gauge link configurations we measure a correlation matrix of 
Wilson loops Ci m (t) for fixed r by smearing the space-like links using numbers n\ 
and n m of spatial HYP smearing iterations. Spatial HYP smearing means that 
only staples restricted to spatial directions are used and we therefore need only two 
parameters, which we set to a 2 = 0.6 and a 3 = 0.3. This second smearing step 
corresponds in the Hamiltonian formalism to the construction of a variational 
basis of operators 0\ that create a state consisting of a static quark and anti- 
quark pair. We use a basis of the operators, labelled by / = 1, 2, 3. The numbers 
of smearing iterations at each level n\ are listed in Table 2. They are chosen 
such that the physical extensions of the operators are approximately constant as 
the lattice spacing changes. Finally we use the generalized eigenvalue method to 
extract the static potential from the correlation matrix C[ m (t), the details of this 
can be found in [36]. 

The solution of eq. (4.1) is found by interpolation of the force F, using a 
2-point interpolation F(r) = fo + /bA" 2 - I n order to control the systematic error 
we compare the result with a 3-point interpolation adding a /4/r 4 term. We find 
the systematic error to be negligible. The error analysis takes into account the 
coupling to the slow modes as explained in Sect. 2.2 and we neglect the systematic 
error due to excited state contributions to the potential, because we ensure that 



r 2 F(r)\ =1.65, 

v / 1 r=ro j 



(4.1) 
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p 


S 2 = 0, c = 1.1 


S 2 = 0, c = 1.4 


c= 1.4 


c= 1.8 


5.2 


6.145(59) 


6.119(45) 


6.22(16) 


6.160(79) 


5.3 


7.259(66) 


7.226(47) 


7.33(17) 


7.274(83) 


5.5 


10.00(11) 


9.953(90) 


10.10(26) 


10.04(13) 



Table 3: Chiral extrapolations of Rq from global fits of the form eq. (4.5). The cut c indicates 
the upper bound in x; the fits with S2 = do not include a quadratic term in x. The final 
results are taken from the first set. 

it is much smaller than the statistical one. The values of Rq for each ensemble are 
listed in Table 1. 

In order to define the lattice spacing in a mass independent way we need to 
perform a chiral extrapolation of R at fixed (5. The chiral fits are done in the 
variable 

x = {r m PS ) 2 . (4.2) 

Before discussing the chiral fits we study the cut-off effects in the mass dependence 
of r . For this purpose we define a reference value r 0rc f, which corresponds to the 
value of r at the pseudoscalar mass 

^ref = (r mp S ) 2 \ rc{ = 0.75 . (4.3) 

The reference point eq. (4.3) corresponds to a pseudoscalar mass of 340 MeV. In 
the left plot of Fig. 3, we plot ro/Vo rc f versus x at our three different (3 values. 
The value of r 0rc f is obtained by linear interpolation and the error analysis of 
W r 0rcf takes into account the correlations between the data. We only consider 
data with x < 1.4. Using these data we determine the first coefficient in the 
Taylor expansion of r /r 0rc f(x) around x TC f independently for each value of j3 

-^-(x) = 1 + s(a/r 0rcf ) ■ (x - x Tei ) . (4.4) 

^Oref 

We do not find significant cut-off effects in the slope s(a/r 0ici ) [39] and in the 
continuum limit we obtain the value s(0) = —0.067(10) by fitting to a constant. 

Motivated by the results in the left plot of Fig. 3, we perform a global fit for 
all (5 values simultaneously of the form 

Ro{fi,x) =R 0c (P) (1 + S 1 x + S 2 x 2 ) (4.5) 

and with cuts 

x < c (4.6) 



10 




Figure 3: The left plot shows the ratio ro/r"o re f as a function of x = (romps) 2 , where roref is 
defined at the reference point eq. (4.3), and the cut x < 1.4 is applied. The right plot shows our 
data of Rq = r^/a with the chirally extrapolated values using a linear function in x (solid lines) 
applying the cut x < 1.1. For comparison a quadratic fit with c = 1.8 is also shown (dashed 
lines) . 

on the pseudoscalar mass. The cuts c = 1.1 , 1.8 correspond to mps = 410 MeV, 
530 MeV. The fit takes into account both errors on R$ and on x. The chirally 
extrapolated values Ro c (f3) are listed in Table 3 for various possibilities of the 
global fit. Columns two and three are linear fits while columns four and five are 
quadratic fits. We quote as our final numbers for R Qc ((3) the results in the second 
column of Table 3 from the linear fit with c = 1.1 (which has x 2 /dof = 0.07 and 
S 1 = -0.073(14)): 

#00(5.2) = 6.145(59) , R 0c (5.3) = 7.259(66) , R 0c {5.5) = 10.00(11) . (4.7) 

The data for Rq(/3,x) and the linear fit with c = 1.1 (solid lines) are shown 
in the right plot of Fig. 3. The numbers in eq. (4.7) cover all the fit results of 
Table 3 within errors. In particular they are perfectly consistent with the results 
in the fifth column of Table 3 from a quadratic fit applying the cut c = 1.8 (which 
has x 2 /dof = 0.28, Si = -0.085(25) and S 2 = 0.013(11)). This quadratic fit is 
represented by dashed lines in the right plot of Fig. 3. The quadratic term is not 
significant even with c = 1.8. 

The discrepancy with the determination of r by QCDSF [40] was discussed 
in [41]. Meanwhile QCDSF updated their values, 4 which now agree with our 
determination where they can be compared. 

4 G. Bali, private communication. 
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5 Chiral extrapolation of /k and strange quark mass 

In this section we describe our central determination of the scale. For a number 
of reasons, it is based on the kaon decay constant /k- First, /k is experimentally 
accessible once the CKM-matrix element V^ s is considered known, which is a good 
assumption within the envisaged precision. Second, chiral perturbation theory 
(ChPT) provides a theory for the quark mass dependence of f K at small masses 
of the light quarks, i.e., our sea quarks. Third, as already mentioned, we remain 
within the pseudoscalar sector of the theory, where ground state properties can 
be determined without doubt. 

5.1 Our strategies 

The main difficulty and source of a systematic error is the extrapolation to the 
proper quark masses, the "physical point". Once we decide to set the scale through 
/k, this point is naturally defined by 

R K = R^ S , R, = Rt y \ (5.1) 

where 

JK JK 

and i?P hys , Rfy* are the values of these ratios in Nature. In an attempt to 
minimize uncertainties, we take the physical masses and decay constants to be 
the ones in the isospin symmetric limit with QED effects removed as discussed 
in [42]. We use 

™* lP bys = 134.8 MeV , m K , P hys = 494.2 MeV , / K , P h y s = 155 MeV . (5.3) 

The two conditions eq. (5.1) define a point in the plane spanned by mnght, m s , 
or equivalently k±, k 3 . We are presently not able to simulate at or very close to this 
physical point, especially not for the smallest lattice spacing, where huge lattices 
would be needed to keep finite size effects under control. The physical point has 
to be approached from unphysically large values of m^t- Along which trajectory 
in the plane of bare parameters one approaches the physical point is in principle 
arbitrary. However, one would like the quantity that is to be computed — here /k 
- to depend very little on the distance to the physical point, allowing for an easy 
extrapolation. Secondly, one would like a theoretically motivated extrapolation 
formula. In ChPT both TOk and /k depend on the sum of the quark masses at 
leading order in the systematic expansion in small quark masses. Keeping 

i2 K («i,« 3 ) = J R? yB (5.4) 
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0.25 
0.2 



0.15 
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Figure 4: Left: Trajectories to approach the physical point in the plane of light and strange 
quark mass -Mn g ht, M a . The dotted line corresponds to strategy 1, i.e., i?x = R^ yS ; whereas 
strategy 2 (full line) holds M s fixed. Right: The two functions Lk and L w , eqs. (5.9) and (5.13) 
respectively, in the interval yi g [0,j/k]- 

thus defines a trajectory where /k varies little in ChPT. We will discuss this quan- 
titatively below. An additional advantage is that all along this trajectory we have 
ffiK ~ ^x hyS ' wnne for a more conventional trajectory, where m s is kept constant, 
the mass itik is significantly heavier than m^ hys . Since the ChPT expansion is 
written in terms of m|- and m^, having mx no larger then m^ hys increases our 
chance of being inside the expansion's domain of applicability. To our knowledge 
this strategy has not been used so far, somewhat surprisingly. 

As an alternative we use a second strategy, where we keep m s constant, 
however instead of using the expansion in and m^, we use SU(2) ChPT for 
kaon observables [43,44], where only is considered a small parameter (in units 
of the chiral scale 47r/). This strategy provides a suitable definition for future 
investigations of mesons and baryons with strange quark content. 

The trajectories belonging to the two strategies, whose details are the subject 
of the following sections, are schematically shown in Fig. 4 on the left. 

Since the following section deals with the scale setting, we want to distinguish 
between the kaon decay constant in physical units /k and in lattice units Fk 

F k = af K . (5.5) 



13 



5.2 Strategy 1 and partially quenched SU(3) ChPT 

As discussed above, the trajectory k 3 = h(Ki) is defined by 

R K ( Kl ,h( Kl )) = 



(5.6) 



where Rk is considered a function R k (k,i,k 3 ) . In practice we take a fixed sea 
quark hopping parameter K\ and a few values of k 3 for which R k (ki,k 3 ) is close 
to i?^ hys and then interpolate in k 3 to find h(ni). Also /k(^i, ^3) is interpolated to 
K3 = H(ki). We give more details on the various interpolations in K3 in Sect. 5.4. 

It remains to extrapolate /k( k 1) m K i to the physical point eq. (5.1). 

We use the partially quenched chiral perturbation theory results by Sharpe [45], 
implement our condition eq. (5.4), which expresses uik in terms of m n , and find 



/k(k1,Mki)) — /K.phys 

Lk(Vi, Vk) = L K (y 1 ,y K ) - L K (y w , y K ) , 



1 + L K (y u y K ) + («4 - |) (yi - y^) + 0(y 2 ) 



Lk{Vi,Vk) = -\vi log(yi) - ^yi log(2y K /yi - 1) 



(5.7) 
(5.8) 
(5.9) 



The variables y, are proportional to (averages of) quark masses up to quadratic 
terms 



Vi 



w4(ki) 



Stt 2 / 2 ^,^))' 



(5.10) 



m 



7r,phys 



87T 2 /, 2 



0.00958 



m 



Vk 



K,phys 



K,phys 



0.12875 



,phys 



Because of eq. (5.4), we have y 3 = m^Ki, /i(/ci))/[87T 2 /k(/Ci, = 2yK — yi + 

0(y 2 ) and y 3 does not appear in eq. (5.7). 

At this order in the chiral expansion we can also replace 

5,1 = iw&rv (5 ' n) 

with the corresponding replacement y^ — > y^, which we will use as a check of the 
typical size of 0(y 2 ) effects. 

In the right panel of Fig. 4 we compare the chiral log function Lk to the one 
describing the chiral behavior of /„., 



1 + L w (j/i) + (a 4 + -a 5 ) (yi - y^) + 0(y 2 ) 



,phys 

L v (yi) = -yi log(yi) , ^(yi) = - L^y*) ■ 



(5.12) 
(5.13) 
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Our condition eq. (5.6) leads to the specific combination of chiral logarithms 
eq. (5.9), which has very little curvature and is overall much smaller than L n ; the 
suppression of the light quark mass dependence thus extends also to the NLO 
chiral logarithms. This suggests that the chiral extrapolation is much easier than 
for /„. and was one of our reasons to select /k to set the scale. Of course, the 
counter terms ctj do not contribute in Fig. 4 (right), but as they are linear in yi 
introduce no curvature. 



5.3 Strategy 2 and SU(2) ChPT 

Here we extrapolate in the light quark mass at fixed mass of the strange quark, 
namely we tune for each sea quark mass the strange quark's hopping parameter 
such that the PCAC mass am 3 4 has a prescribed value /i, which is independent of 
K\. b This defines the function k 3 = s(ki,/i). In practice, we again interpolate the 
data for m 34 (/ti, k 3 ) in k 3 and then solve 

aTO 34 (/ci,s(«i,/i)) = /i (5-14) 

for s, with the left hand side represented by the interpolation formula. 

To find the value of /i corresponding to the physical point, we employ SU(2) 
ChPT [43,44] to first extrapolate M| = (am K ) 2 and F K (both interpolated in k 3 
to the point k 3 = s(ki, //)) in y\ to y\ = y^ at fixed value of //, 

3 

F K («i,s(«i,ji)) =p(ji) [l - g[yilog(yi) -y*\og(y w )] +a t (fi) (yi - y v ) + 0(yf)] , 

(5.15) 

Mk(ki, s(ki, fi)) = q(fi) [1 + a m (fi) (yi - y*) + 0(yf)] . 

These expressions represent the asymptotic expansions for small y\ at fixed \x 
correct up to error terms of order yf. 6 

From eq. (5.15), q(n) and p(fj) are computable for arbitrary values of \x. The 
requirement that ^k//k attains its physical value at the physical light quark mass 
then defines /i s , 



2 



^£ (5.16) 



5 In principle we should keep m^ 4 fixed, but the ratio m^/m^ is independent of Ki since we 
set 6a = bp — 0, see eq. (3.4). 

6 We note that SU(2) ChPT for kaons does not take into account kaon loops as opposed to 
SU(3) ChPT. This corresponds to the production of two kaons, i.e., states with an energy of 
around one GeV. From this point of view it is not really worse than ChPT in the pion sector 
where \pw) intermediate states are dropped. We thank Gilberto Colangelo for emphasizing this 
point. 
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This equation is solved numerically for /i s and the lattice spacing is then given by 

Pith) 
a = - 

JK.phys 

As before, the constants ccf(// s ), a m (/x s ) have common values for all three (3 in the 
fits to eq. (5.15). 

5.4 Interpolations 

In various places we need observables, such as Rk(ki,K3) as a continuous func- 
tion of k 3 , not just for a few numerical values. In all cases we chose four dif- 
ferent numerical values for n 3 , close to the required one, namely those which 
have the smallest distance defined by d — \Rk(ki,K3) — R^ ys \- We then de- 
termine an interpolation polynomial P2(«3 — «i) = Yln=o Cn ( K 3 ~~ Kl ) n through 
a fit to the data with weights (d 2 + e 2 ) -1 . The regulating term e is chosen as 
e = (max K3 {i? K ( K i ? %)} ~ mm K 3 {FK(/ti, k 3 )}) • We always checked that un- 
certainties due to the specific interpolation are negligible by considering natural 
variants. 

5.5 Cutoff effects 

So far the discussion neglects lattice artifacts completely. We will see that these are 
very small in our formulation. Since both y\ and a 2 are small, it is natural to keep 
terms of order y 1 (as done above) and those of order a 2 which are independent of y\ 
but drop terms of order y 2 , yia 2 , a 3 etc. In this approximation we can use eq. (5.7) 
for the decay constants in lattice units Fk, and with a global, a-independent low 
energy constant a 4 . The straight a 2 term does not contribute to eq. (5.7), since 
/x.phys is used to set the scale. Such a term is, however, present in general. An 
example is the combination Tq/k- It is shown in Fig. 5 with open symbols, where 
both / K and r are evaluated at the finite quark mass. Since a ChPT expansion 
of r does not exist, we use 

i2o(«i) F k (ki, = [ro\y=v* /K, P h y s] cont [l + A ( Vl - y w ) + BF 2 ] (5.18) 

to describe the data and to extrapolate to the physical point and continuum limit 
[ro\ y=y7I /k,ph y s] cont - Taking into account that we have < 4 • 1CT 3 , the lattice 
artifact parametrized by B is rather small, see Table 4. In these fits we take the 
chirally extrapolated value Fk in the last term of eq. (5.18). Replacing it with 
Fk at the given value of yi, we get an insignificant shift of 0.6% in the continuum 
value, which tests the smallness of yia 2 terms. 
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(5.17) 



„,max 
V\ 


r o|2/=j/„-/K,phys 


A B 


?"0c./k,phys «4 


B' 


0.10 


0.3951(62) 


0.28(22) -20(5) 


0.3950(75) 0.52(11) 


-22(6) 




Table 4: 


Results for chiral fits in 


the range y± < j/™ ax . 





0.5 



0.45 - 



o 



0.4 



0.35 



□ f3 = 5.5 
o 13 = 5.3 
A 13 = 5.2 




Vtt 



0.15 



Figure 5: Chiral extrapolation of tq/k- For the open symbols ro and /k are evaluated at finite 
quark mass, the dashed lines show eq. (5.18) for the three values of (3. The filled symbols use 
the extrapolated value r = r 0c , the dashed dotted lines represent eq. (5.19). In both cases, the 
corresponding solid line gives the continuum result. 



Alternatively, we use the already chirally extrapolated values Rq c from Sect. 4 



and 



#oc-Fk(/Si,M/si)) 

-physj cont 



1 + L K (yi,y K ) + (a 4 - j) (j/i - y v ) + 



. (5.19) 



The dashed dotted lines in Fig. 5 show the corresponding fit with the parameters 
listed in Table 4. 

The just discussed fits mainly illustrate that cut-off effects are small and 
provide a motivation that indeed terms such as a 2 y\ can be dropped. We then 
determine the values of the kaon decay constant in lattice units at the physical 
point .Fk.phys by a direct application of 



Fk(«1, h{Kl)) — -Fk.phys 



1 + L K (yx,y K ) + (o4 - -) {y x - y n ) 



(5.20) 



Here the data at the three different f3 are combined in a global fit with one common 
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value of o?4, but, of course, with different F Kphys for the three different f3. The 
lattice spacings are then obtained from 



a= F p^. (5.21) 

jK.phys 

Results are discussed in Sect. 5.6 together with our alternative strategy. 

We finish the discussion of cut-off effects by mentioning an investigation of 
the PCAC relation. In a properly 0(a) improved theory we have A = \{rn^ + 
m^ 4 ) — m^ 3 = 0(a 2 ). We have computed A for all our combinations of (3, k±, k 3 
and found aA < 1CT 4 and A < 0.3 MeV in physical units. 7 This confirms again 
the smallness of cut-off effects. 



5.6 Numerical results 

We now apply the above formulae to a determination of the kaon decay con- 
stant and the strange quark mass. The renormalization of the decay constant, 
eq. (3.8), starts from Za computed in [32] and improved in statistical precision 
in App. B. For the mass-dependent 0(a) improvement terms, which yield very 
small corrections, we use 6a = and 6a = 1 + 0.06167gg from one-loop pertur- 
bation theory [31]. The extrapolations of F K through eq. (5.20), strategy 1, and 
eq. (5.15), strategy 2, are shown in Fig. 6. In the case of strategy 2, the parameter 
/i is fixed to the strange quark mass /i — fj, s through eq. (5.16), with the values 
of k s given in Table 1. The fits include data for yi < y™ ax with y™ ax = 0.1. We 
observe that the two rather different chiral effective theory extrapolations yield 
results in close agreement at the physical point. Table 5 lists the fit parameters 
and also includes fits with y\ — > y±, simple linear fits in y\ and y\ as well as results 
including data out to y™ ax = 0.15. 

A few observations are worth pointing out. Using y\ as the chiral variable 
leads to a stronger dependence on the cut ?/™ ax . Still, the extrapolation in terms of 
yi shows a tendency to converge to the one in terms of y\ when the cut is lowered, 
as it should be. Results of linear extrapolations are also very close. 

We therefore take our central values from the results of strategy 1 with cut 
^max = o.l. As a systematic error we take into account the deviations to the fit 
following strategy 2 and to a simple linear extrapolation; for our smallest lattice 
spacing these different extrapolations are compared more closely in Fig. 7. The 
final numbers for Fk at the physical point and the associated lattice spacings are 
shown in Table 6. 

7 The improvement term proportional to bj± — bp = (bj^ — bp)/Z was inserted non- 
perturbatively using 6a — 6p and Z from [46]. It is very small in practice. 
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Strategy 2 




10 2 F K 




10 2 a[fm] 


10 2 F K 


at 


10 2 a[fm] 


(5.15) 0.15 


5.2 


5.86(7) 


1.30(6) 


7.45(9) 


5.77(7) 


1.42(6) 


7.35(9) 




5.3 


5.09(5) 




6.48(7) 


5.01(5) 




6.38(7) 




5.5 


3.76(3) 




4.79(4) 


3.71(3) 




4.72(4) 


0.1 


5.2 


5.88(7) 


1.13(8) 


7.49(9) 


5.83(7) 


1.15(8) 


7.42(9) 




5.3 


5.11(5) 




6.51(7) 


5.06(6) 




6.44(7) 




5.5 


3.79(3) 




4.82(4) 


3.75(3) 




4.77(4) 


linear 0.1 


5.2 


5.97(7) 


1.78(8) 


7.61(9) 


5.92(7) 


1.72(8) 


7.54(9) 




5.3 


5.19(5) 




6.61(7) 


5.15(5) 




6.55(7) 




5.5 


3.84(3) 




4.89(4) 


3.81(3) 




4.85(4) 



Table 5: Results for chiral fits in the range y 1 < y" 



The fit parameters and «f turn into low energy constants of SU(3) and 
SU(2) chiral perturbation theory, respectively, when the cuts on y are sufficiently 
small, such that the used NLO expressions are accurate. A look in the tables 
shows that their dependence on the cut is anything but negligible. Clearly the 
values obtained with y™ ax = 0.1 are our best estimates. Differences to y™ ax = 0.15 
are an estimate of the systematic error, which, however, deserves a more detailed 
study. We postpone this to a later analysis when data at an even smaller quark 
mass will be available for two values of the lattice spacing. 

6 Lambda parameter 

The A parameter, which normalizes the scales in the perturbative running of QCD 
quantities at high energies, has been determined in an earlier work by the ALPHA 
collaboration [5] in terms of a length scale L 1 . It used the Schrodinger functional 
which allows to compute the running over many orders of magnitude in the scale 
L, and L\ is a typical hadronic scale which allows to make contact to physical 
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p 


F K 


a[fm] 


5.2 


0.0593(7)(6) 


0.0755(9)(7) 


5.3 


0.0517(6)(6) 


0.0658(7)(7) 


5.5 


0.0382(4)(3) 


0.0486(4) (5) 



Table 6: The kaon decay constant at the physical quark mass and the corresponding lattice 
spacing. The first error is statistical, the second systematic error due to the chiral extrapolation. 
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Figure 6: Chiral extrapolation of the kaon decay constant in lattice units for all three j3. Open 
symbols and dashed lines correspond to strategy 1, whereas filled symbols and dash-dotted lines 
represent strategy 2. 



units. 

In particular a renormalized coupling g(L) can be defined for which the renor- 
malization group invariant A parameter is given by 

LA = (bofiL^/^e^/^^exp 

Here /3(g) = —Ldg(L)/dL with bo and b\ the first two universal coefficients of 
its asymptotic expansion /3(g) = —g 3 (bo + b{g 2 + ...). Then A is independent of 
the scale but dependent on the renormalization scheme adopted for the coupling. 



9{L) 



dx 



(3(x) ~^ b x 



blx 



■(6.1) 
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Figure 7: Detailed comparison of strategy 1, strategy 2 and a linear extrapolation at = 5.5. 
The symbols are the same as in Fig. 6 with the solid line representing the linear extrapolation 
of the open symbols. Note that the slopes given by the low energy parameters a are determined 
from the global fit to all three j3. 

The relation of the Schrodinger functional A to its value in other schemes is given 
exactly by a one-loop perturbation theory computation. 

For the scale at which we make contact to physical units we use L\ defined 
through g 2 (Li) = 4.484. At this value of the coupling we have a large and precisely 
tuned set of pairs (Li/a, 0). Reanalyzing the data of Ref. [5] at this point, we get 
the continuum value of ALi = 0.264(15). What remains to be done is to compute 
JkL\ in the continuum limit. To this end, we have to combine the updated data 
for Li/a of App. D with the decay constant at physical quark masses described 
in the previous section. Since the data points for L\ja and are at different 
values of /3, we need to interpolate L\ja to the values of the latter; then the 
product /k-^i can be extrapolated to the continuum limit. 

6.1 Interpolation of L\ 

We start with values of L x /a at (3 = 5.2, 5.2638, 5.4689, 5.619, 5.758, 5.9631 
obtained in App. D and Ref. [47]. Since L\ is a physical scale, we expect log(Li/a) 
to be roughly linear in /3. Four different fits have been tried: a linear fit to the full 
range in (3 and to the range [5.2,5.619], which we compared to quadratic fits in the 
same intervals. The central values derived from these interpolations differ only by 
a small fraction of the assigned statistical uncertainty. Also the statistical errors of 
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p 


Li/a 






m s /f K 


5.2 


5.367(82) 


0.318(6)(3) 


1.155(22) 


0.530(12)(6) 


5.3 


6.195(51) 


0.320(5)(4) 


1.169(15) 


0.577(H)(7) 


5.5 


8.280(80) 


0.316(4)(2) 


1.213(17) 


0.617(H)(5) 


cont. 




0.315(8)(2) 


1.252(33) 


0.678(12)(5) 



Table 7: Values of L\/a, Li/k, tq/L\ and to s //k at the three values of /?. For the latter three, 
we also give the value extrapolated to the continuum limit. The running mass in the Schrodinger 
Functional scheme fn s is given at the renormalization scale L\. Statistical and systematic errors 
are given. 

the interpolated points are virtually equal for these procedures, only the quadratic 
function to the restricted interval gives larger uncertainties, as expected from a 
fit with three parameters to four points. Since we want to obtain conservative 
estimates of the errors, we take the result of the latter in the further analysis. 
The three values of Li/ a which we use in the following are listed in Table 7. 

6.2 Continuum value of the A parameter 

Fig. 8 shows that in the continuum extrapolation of Li/k the cut-off effects are 
smaller than the statistical uncertainties; the numerical values are given in Table 7. 
Indeed, a constant extrapolation gives a xVd-o.f. below unity. However, we 
use a linear extrapolation, allowing for 0(a 2 ) effects hidden by the statistical 
fluctuations. 8 In this extrapolation the covariance between the three points from 
the interpolation in Li/ a is taken into account. For comparison with results in 
the literature, we give all results also in units of r \ yi=y ^. The two combinations 
evaluate to 

f K Li = 0.315(8)(2) r /Li = 1.252(33) (6.2) 

in the continuum limit. 

As in the previous section, the final results come from strategy 1 for the 
chiral extrapolation of Fk, given in Table 6. Strategy 2 is used to estimate the 
systematic uncertainty; however, it is small compared to the statistical errors and 

8 That the leading corrections come at 0(a 2 ) is based on the assumption that the two-loop 
approximation of the boundary improvement coefficient ct of the Schrodinger functional is suf- 
ficient at the accuracy we are aiming at. This is corroborated by experience from the quenched 
approximation [48] and earlier investigations of the approach to the continuum with c t at 1-loop 
and 2-loop precision [49]. 
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Figure 8: Continuum extrapolation of L\ in units of /k. Even though the data shows no cut- 
off effects, we use a linear extrapolation to account for uncertainties from 0(a 2 ) effects hidden 
by the errors. The two strategies for the chiral extrapolation of Fk = a/x agree well within 
statistics. 

we do not give it separately. We therefore quote 



Now, as a result of our analysis, the error is dominated by the error on AL\. We 
translate to the MS scheme using AH = 2.382035(3)A^ [2,50] and find 



7 Strange quark mass 

The determination of the strange quark mass uses a strategy [51] which splits up 
the computation into several steps, circumventing the multi-scale problem a direct 
approach would face. It therefore allows for a good control over the systematic 
errors. Earlier results have been presented in Ref. [6], here we improve on them 
with a much higher statistical and systematic accuracy. 

To briefly summarize the procedure, the RGI mass 9 and the renormalized 
running mass m s in the Schrodinger functional scheme are given in terms of the 

9 The RGI mass can be defined through M s = lim L ^ (2& g 2 (£)r do/(2bo) m s (X) with the 
universal coefficient do = <V(47r) 2 [52]. 



A( 2 )// K = 0.84(6) and r A {2) = 0.331(22) . 



(6.3) 



A|| = 310(20) MeV and r A^ = 0.789(52) . 



(6.4) 



23 



bare PC AC mass m s by 

M s = Z M m s = —-m s {L) = —-^—-m s , 7.1 
m(L) m(L) Zp(L) 

where the computation of the renormalization constant Zm is split in two steps. 
The continuum value of the universal first factor M/fn has been computed in 
Ref. [6] for the two flavor theory which we consider here. It does not depend 
on the quark flavor or the lattice regularization. For the second factor, updated 
values of Z P are presented in App. C. Of course, the first factor and Z P have to 
be evaluated at a common renormalization scale L. While Ref. [6] used L = L max 
with g 2 (L max ) = 4.61, we choose here L — L\ with g 2 (Li) = 4.484 as in the 
previous section. The first factor in eq. (7.1) is M/fn = 1.308(16) at this scale 
from a re-evaluation of the data of Ref. [6]. 

7.1 Determination of the strange quark mass 

For the bare strange quark mass at physical light quark masses m s , the final 
ingredient of eq. (7.1), we again follow the two strategies presented in Sec. 5. 
Of the two, the second holds the strange quark mass fixed as a function of the 
light quark mass and is therefore the natural choice for this analysis. We still 
use strategy 1 as a cross check to estimate systematic effects from the chiral 
extrapolation. 

7.1.1 Strange mass from strategy 2 

The hopping parameter of the strange quark k s is a non-trivial function of Ki, 
defined in strategy 2 as 

« s = s(ki,^b), (7.2) 

with s(ki, [/,) defined in eq. (5.14). Its numerical determination as well as the 
determination of (j, s have been explained in Sect. 5.3. The resulting values of k s 
are listed in Table 1. This definition holds the bare PCAC strange quark mass 
fixed to the value /i s set by the physics requirement eq. (5.16) at physical light 
quarks. 

Expressing M s in units of /k we eliminate Z\ and get 

M k^[1 + (>a->pM^, (7.3) 



/K, P h„ m(L) Z P (L) 1 ^ r,m F*T 



K,phys 



with -^xfphys = ^k,phys/^A- The second factor is 0(a) improved, if we neglect, as 
before, a tiny correction (6a — bp)am sea . 
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7.1.2 Strategy 1 

Within this strategy, it is most natural to determine the combination M s + M, the 
sum of strange and light quark mass, since this combination is kept fixed at lowest 
order in chiral perturbation theory along the trajectory defined by eq. (5.6). The 
first order corrections are easily incorporated from [45] such that we arrive after 
renormalization as in eq. (7.3) at 



2m l3 (K 1 ,h(K 1 )) _m s + m 

,phys 



1 + L m (y 1 ,y K ) + (a 4>6 - -) {y 1 - y n ) + 0(y 2 ) 



(7.4) 

L m (yi,yK) =L m {yi,yK) - L m (y n ,y K ) , «4,6 = 3a 4 - 4a 6 , (7.5) 

3 

L m (yi, y K ) = -(vk- gS/i) log(2y K /yi - 1) - y K log(yi) , (7.6) 

which we can use to extract (M s + M)// K , P hys and the combination of low energy 
constants a 4 6 . Removing the contribution of the average light quark mass M 
amounts to multiplying by a correction factor such that M s = (M s + M)(l — p), 
where p is a small number and can therefore be incorporated by its lowest order 
chiral perturbation theory estimate 



M m 



M s + M 2mi 



= 0.037. (7.7) 



Because of this last approximation, we prefer the determination of M s from strat- 
egy 2 and use the one here just as a consistency check. 



7.2 The strange quark mass in the continuum limit 

For strategy 2, the strange quark mass is the strange quark mass is directly kept 
fixed at its physical value. The discussion can be found in Sec. 5. For strategy 1, 
the extrapolation follows the same principles as before: we neglect cut-off effects 
in the NLO terms of the chiral expansion. The corresponding fit to the data is 
displayed in Fig. 9 on the left, which demonstrates that this assumption holds well 
within the statistical accuracy. Also shown for comparison is a linear extrapolation 
which agrees within the uncertainties with the values obtained from the ChPT 
formulae. 

The values of m s / f K for the two strategies as a function of the lattice spacings 
are plotted in Fig. 9 on the right. Already at finite lattice spacing the two sets 
seem to be shifted with respect to another and therefore have parallel continuum 
extrapolations which, however, agree within uncertainties. The data fits well 
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Figure 9: Left: Chiral extrapolation of the strange quark mass according to strategy 1. The 
dashed line indicates a corresponding linear fit which gives a value that agrees within uncertain- 
ties at the physical point. Right: Continuum extrapolation of the strange quark mass in units 
of /k for the results of the two strategies. 



the assumption that only leading cut-off effects are present and we observe a 
correction of 8% from our value on the finest lattices. Since we prefer strategy 2 
for conceptual reasons, we quote as a final result with statistical and systematic 
uncertainties 

m B // K = 0.678(12) (5), M s /f K = 0.887(19) (7) , M s = 138(3) (1) MeV, (7.8) 

where we use M/fn = 1.308(16) as quoted before. For reference, we also give the 
numbers in the MS scheme. This conversion is the only part of the computation in 
which we need to take recourse to perturbation theory, known in this case to four 
loops, which differs from the two- and three- loop result by only a small amount. 
We use the same method as described in Ref . [6] , but with the new value of Aj^g 
which leads us to m m (2GeV)/M = 0.740(12) and 

mf (2GeV) = ^ m^(2 GeV) = 102(3)(1)MeV (7g) 

Here also the statistical uncertainty of A^jg is included. 
8 Conclusions 

Setting the scale in lattice computations requires a dimensionful quantity which 
is known from experiment to good accuracy for which we chose the kaon decay 
constant. For currently available lattices this involves chiral extrapolations to the 
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physical light quark masses, which can be important corrections depending on the 
quantity. Chiral perturbation theory turns out to be a good guide, if one stays 
within its domain of applicability. In particular the first strategy discussed in 
this paper, designed to have small chiral corrections for the kaon decay constant, 
has proven to be successful in the sense that the extrapolation is very robust 
under reasonable changes in the functional form: various next-to-leading order 
expansions and also a linear extrapolation give compatible results. 

Whereas strategy 1 keeps the sum of strange and light quark mass approx- 
imately constant, strategy 2 provides a constant strange quark mass, making it 
suitable for the determination of M s and further analysis of strange baryons. This 
flexibility is probably particular to the situation of a quenched strange quark, 
however, the fact that the results from both strategies agree within small uncer- 
tainties gives rise to the expectation that also with a dynamical strange quark and 
the u/d quarks sufficiently in the chiral regime, the extrapolations can be safe. 

As already mentioned in the introduction, this publication completes the 
iVf = 2 part of a long research program by the ALPHA collaboration, which has 
the determination of the A parameter and the strange quark mass as two of its 
main goals. We can therefore quote these two observables for the two flavor theory 
in physical units 

A^ = 310(20) MeV, M s = 138(3) (1) MeV , (8.1) 

with the value of the strange quark mass in the MS scheme given in eq. (7.9). 

Since previous results in physical units were wrong due to inaccurate values 
for r Q /a in the literature, we have also computed this scale and get 

r = 0.503(10) fm (8.2) 

at the physical pion mass using f K = 155 MeV. 

The restriction to two flavor QCD introduces some ambiguity to the scale 
setting and quark mass determination, because mass ratios do not have to be the 
same as in Nature. The effect of the heavier quarks is unknown and interesting in 
its own right. Comparing to the Nf — 2 results presented here, future calculations 
will shed light on the question of the contributions of strange and charm quarks 
in the sea. 

For now we can make a comparison to the values obtained in the quenched 
approximation. In the strange quark mass we find no significant difference to the 
corresponding value quoted in Ref. [54]. However, the A parameter exhibits an 
interesting dependence on Nf, with the five flavor theory A^ 5 ^ = 212(12) MeV ex- 
tracted from Bethke's world average a^jg(Mz) = 0.118(10) [53] and the quenched 
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Nf 




experiment 




theory 





238(19) MeV 


m K , K -)• fiu^, K -)• 




lattice gauge theory [3] 


2 


310(20) MeV 


m K , K -)• /xi/^, if ->■ 




this work 


5 


212(12) MeV 


world average 




perturb, theory [53] 



Table 8: The A-parameter as a function of the number of flavors. 

value A(°) = 238(19) MeV [3] listed in Table 8. Perturbation theory predicts 
S_( N t) < A^ -1 ), which is valid and precise when the quark being decoupled is 
sufficiently heavy. This is expected to be the case for the b-quark and thus 
y\( 5 ) < A^ 4 ). However, comparing the ALPHA value for Nf — and our new 
Nf — 2 result shows that for light quarks the A f -dependence is in the opposite 
direction. The A parameter increases quite strongly with increasing Nf. An even- 
tual decrease is likely to be present between Nf — 3 and Nf = 4, leading to an 
agreement with the phenomenological determinations of a^(M z ) [53]. 

In any case, it will be an important milestone to arrive at a non-perturbative 
and precise determination of A^ 4 ^ and compare with the present world average. 
This is particularly needed in the light of the spread of different determinations 
before taking the world average. For example the authors of Ref. [55] analyse deep 
inelastic scattering data by a global fit to parton distribution functions and ay^ 
at NNLO. Their value for Oyg(Mz) corresponds to A^ = 160(11) MeV, which is 
quite a bit smaller than the world average. 
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r 


N 


m 


N 2 


acc. rate 


MD time 


A3 


8 4 


0.37 


2 


4 


5 


125 


0.91 


8030 


A4 
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0.37 
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4 


5 
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0.85 


8090 


E5f 


8 4 


0.37 
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4 


5 
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0.87 


16000 


E5g 
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4 


5 
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0.83 


16180 


F6 


8 4 


0.37 


2 


4 


5 


260 


0.89 


4800 


F7a 


8 4 


0.37 


2 


4 


5 


350 


0.84 


5600 


F7b 


8 4 


0.37 


2 


4 


5 


350 


0.87 


4000 


N4 


8 2 12 2 


0.44 


0.5 


4 


5 


24 


0.88 


3700 


N5 


8 2 12 2 


0.44 


0.5 


4 


5 


24 


0.87 


3800 



Table 9: Parameters of the DD-HMC algorithm for the ensembles used in this publication. We 
give the HMC block size, the corresponding ratio of active links i? ac tivc, the trajectory length, 
the (relative) step sizes of the gauge, block fermion and global fermion force. This is followed 
by the measured acceptance rate and the total statistics after thermalization. 

A Details of the update algorithms 

In our simulation we first used the DD-HMC algorithm [14, 15] as implemented 
in Luscher's package [16] and then switched to our implementation of a mass pre- 
conditioned HMC [23]. Both algorithms have several parameters which influence 
their performance and we therefore give the essential ones in this section. 

A.l DD-HMC 

The block decomposition of the DD-HMC algorithm separates the infrared from 
the ultraviolet modes. The details and notation can be found in Ref. [14]. We give 
the size of the blocks and the step sizes of the three different forces in Table 9. 
The gauge force has NqNiN 2 steps per trajectory, the block fermion force NiN 2 
and the global fermion force N 2 steps. The number of trajectories times their 
length is given in the "MD time" column of the table; to get an effective statistics, 
this can be multiplied by the ratio of active links -R ac tive- 

A.2 MP-HMC 

In our implementation of the mass preconditioned HMC [23], we split the fermion 
determinant detQ 2 , with Q = a^D, first by symmetric even-odd preconditioning 

det Q 2 = det 2 Q cc det 2 Q GO det 2 Q s (A.l) 
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K2 
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N 




N 2 


N 3 


acc. rate 


MD time 


A5c 


0.135887 


0.134250 


2 


9 


i 


1 


32 


0.93 


1160 


A5d 


0.135887 


0.134250 


2 


9 


i 


1 


32 


0.92 


1700 


N6 


0.136552 


0.133857 


2 


9 


i 


1 


16 


0.84 


4000 


07 


0.136550 


0.135000 


2 


9 


i 


2 


16 


0.83 


4000 



Table 10: Parameters of the MP-HMC algorithm for the ensembles used in this publication. 
The Hi parametrize the preconditioning masses and are followed by the trajectory length and 
the (relative) step numbers per trajectory. Also the acceptance rate and the statistics after 
thermalization are given. 

and then apply two levels of mass preconditioning [21] to its Schur complement 

Qs = 1 — QeeQeoQooQoe , 1-e., 

det Ql = det 2 \W~\a 2 )W{a 1 )] det 2 [W _1 (<7i)Qs] det 2 [W(a 2 )\ (A.2) 

with W(a) = Qs+°" an d a 2 > ai > 0. The forces from the individual contributions 
to the action are integrated on multiple time scales [56], the gauge force is on the 
finest level, the determinants of Q ee and Q oa are integrated together with the third 
determinant of eq. (A.2), the next two determinants constitute the two coarsest 
levels. For each we use a second order integrator [57], where an elementary step 
of size e is given by 

T(e) = T 1 (Ae)T 2 (e/2)Ti((l-2A)e)T 2 (e/2)T 1 (Ae) . (A.3) 

T\ is the force application on each level, whereas T 2 is again given by such a step 
with the force of the lower level. We use A = 0.19 for all levels. Only the gauge 
force is integrated with a standard leapfrog. Analogously to the DD-HMC, we 
give the relative step sizes in Table 10, note however, that eq. (A.3) implies a 
factor of two in the number of force applications between levels even for iVj = 1. 
The preconditioning masses are given in terms of the hopping parameter such that 
2a i = re, 1 — k sc \. 

B Improved determination of Za 

The renormalization factor Za was previously computed [32] from a carefully 
chosen Ward identity. At p — 5.2 an error of 2.1% resulted from the fact that at 
this (3 an extrapolation in the quark mass was necessary. Since this uncertainty 
directly propagates into the strange quark mass determination and scale setting, 
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L/a 


T/a 


K 


Za 


am 


5.20 


12 


18 


0.13550 
0.13560 
0.13570 
0.13580 


0.7853(15) 
0.7792(20) 
0.7767(20) 
0.7757(27) 


0.02082(20) 
0.01700(20) 
0.01345(22) 
0.00892(15) 










0.771(6) 


-> 



Table 11: Simulation parameters and results for Za and the PCAC mass am at /? = 5.2. The 
last row shows the chirally extrapolated value as described in the main text. To be compared 
with Table 1 of [32] and its chiral value 0.774(16). 

we here improve Z A at j3 = 5.2. We substantially increased the statistics for the 
three previously simulated quark mass points and added another mass point, in 
order to better control the subsequent chiral extrapolation. We follow the lines 
of [32], i.e., in the determination of Z A we use the optimal wave function of [29] 
to suppress the first excited state contribution in the pseudoscalar channel, use 
the "massive" definition of the renormalization constant [58], and we drop the 
disconnected quark diagrams which only amount to 0(a 2 ) effects on Z A . In the 
notation of [32] we thus compute Z A on with u = co n (o). The results are listed in 
Table 11. 

We have performed several chiral extrapolations, including weighted averages 
of the two and three data points with the smallest masses, but in the end all 
resulting central values at am = are covered by our preferred ansatz, a linear 
extrapolation of the data points with am < 0.02. It is shown in Fig. 10. The 
total range of the y-axis of the figure is the previous [32] ±lcr error. We now have 
a considerably improved value of Z A = 0.771(6). What remains to be done is to 
provide a global fit for Z A {g^). Again, several fit ansaetze have been tried and all 
of them give comparable results, especially in the /3-range which is important in 
the present paper. Constraining the asymptotic behaviour to the perturbatively 
known one-loop estimate, we finally obtain 

Z A {gl) = 1 - 0.116458 g 2 + c l9 4 + c 2 g$ , (B.l) 

/+1.16\ ift _ 2 ^ . , n ift 4 /+0.74 -0.75\ 

where Cov(c)/3 is the covariance matrix of the parameters c±, c 2 obtained directly 
from the statistical errors of the data points. Due to the constraint fit ansatz 
with just two parameters, the error of the fit is significantly smaller than that of 
the individual data points. We therefore recommend to conservatively use Cov(c) 
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Figure 10: Data points with highly improved statistics for Zx at /3 = 5.2. The point in the 
chiral limit is obtained by a linear extrapolation of the three lightest PCAC masses. 

with the factor 3 applied as above, which yields the grey error band in Fig. 11. 
C Determination of Zp 

From the simulations reported in [47] we also have values for Zp, the renormal- 
ization constant of the pseudoscalar density. They are based on simulations with 
8000 essentially uncorrelated measurements and thus have a very small statistical 
uncertainty. However, the small mismatch of the simulated {L/a, (3, am) has to 
be accounted for and introduces additional errors. We need to distinguish /3 = 5.2 
and (3 > 5.2638, since in the latter case, we only have very small shifts. 

We start with j3 = 5.2. Here we have two different sea quark masses each 
at L/a = 4,6 (see also Table 13). They enclose Lm = 0. We thus first perform 
an interpolation to vanishing mass and subsequently interpolate In Zp to our tar- 
get coupling linearly in Ing 2 , c.f. Ref. [6]. A systematic error is estimated by 
comparing to an interpolation linear in g 2 . 

Now we turn to (3 > 5.2638. Both the shifts to am = 0, g 2 = 4.484 and the 
propagation of the errors in g 2 (L) of Table 13 is done by a first order expansion 
of 

F(x,y) = lnZ P |^ =const with x = ln(<f ), y = Lm , (C.l) 

in x, y. The derivative d x F is computed from perturbation theory, using the 
known (continuum) two- loop and three-loop perturbative results for r(g) and 
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Figure 11: Result for our global fit of Za(3§), together with the error band from eq. (B.l). 
For comparison we also show the perturbative (1-loop) result for Nf — 2. 

0(g), respectively. Note that perturbation theory is known to be rather accurate 
for the running of the mass (and therefore Zp) in the SF-scheme [6]. The other 
derivative is estimated as follows. At (3 = 5.2, L/a = 4, 6 we have two mass points 
each and estimate d y F for g 2 = 3.80, 4.95 by finite differences. These are then 
interpolated to <9 y F(hi(4.484), 0) = 0.148(30). Up to tiny effects coming from 
the renormalization of m, the derivative d y F is universal. We neglect the a 2 /L 2 
corrections to exact universality, using d y F(ln(4A8A), 0) = 0.148(30) at all /3. 

The results for Zp are collected in Table 12. Apart from the special case 
f3 = 5.2, where the interpolation in g 2 introduces a noticeable uncertainty, all 
errors in the last column are dominated by the statistical errors in g 2 (Li). 

Within our uncertainties the ^-dependence of Zp at fixed L = L\ is not 
clearly visible. We thus take the average value of all six results and assign an 
error of about 1% to cover all results: 

Zp = 0.5184(53) , 5.2 < (3 < 6.0 . (C.2) 

D Improved determination of L\ 

The practical reason for g 2 (Li) = 4.484 instead of g 2 (L max ) = 4.61, chosen in 
[5], is that meanwhile dedicated tuning runs were carried out [47], which fix the 
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L/a (3 k am Zp(L,/3,am) Zp(L±,P,0) 



4 


5.2 


0.135041 


-0.00279(14) 


0.57672(15) 


— 






0.134800 


+0.00464(20) 


0.58291(22) 


— 








-> 


0.57904(13) 




6 


5.2 


0.135617 


-0.00747(20) 


0.48796(28) 


— 






0.135226 


+0.00571(17) 


0.49641(24) 










-»• 


0.49275(18) 






5.2 









0.5237(45) 


6 


5.2638 


0.135985 


-0.00585(18) 


0.50676(28) 


0.5150(16) 


8 


5.4689 


0.136700 


-0.00339(13) 


0.51103(29) 


0.5170(17) 


10 


5.6190 


0.136785 


-0.00260(9) 


0.51255(30) 


0.5186(25) 


12 


5.7580 


0.136623 


+0.00040(6) 


0.51856(29) 


0.5187(43) 


16 


5.9631 


0.136422 


-0.00107(4) 


0.51879(35) 


0.5172(48) 



Table 12: Pseudoscalar renormalization constant Zp from simulations in the Schrodinger func- 
tional. The kinematical setup is different from those in Table 13 and thus the PCAC mass 
am differs. The rightmost column lists the results after correcting for an imperfect tuning as 
explained in the text. For /3 = 5.2, two non-zero quark masses have been simulated, allowing 
for an interpolation to vanishing mass. 

simulation parameters (L/a, (3, K sea ) at the line of constant physics 

g 2 (L x ) = 4.484 , L im = . (D.l) 

Here m is the PCAC mass of the degenerate sea quarks and additional details 
can be found in the appendices of [47]. Of course fixing such a line of constant 
physics is only possible within a certain (statistical) uncertainty. While in the 
results of [47] the error due to the line of constant physics is not so relevant, we 
here aim for a more precise setting in order to extract the A parameter and the 
scale L\ in physical units. For this reason we have set up new simulations for g 2 at 
f3 G {5.2, 5.2638, 5.4689, 5.6190}, corresponding to the first five sets of parameters 
in Table 13 which enclose the /3-values used here. For completeness we also list the 
additional two sets for (3 G {5.7580, 5.9631}, taken over from the aforementioned 
reference. The simulations for g 2 proceed along the lines of [5]. Results are listed 
in Table 13 together with the deviation 5[g 2 ] = g 2 (L) -4.484. At (3 = 5.2 there is 
a large deviation since our code does not allow to simulate L/a = 5; we have to 
rely on neighbouring points with L/a = 4,6. At each of these values of L/a we 
have simulations at slightly positive and negative m from which we interpolate to 
Lm = 0, see Table 13. We finally just take the value from L/a = 6 into account, 
since it is closer to our target g 2 . At fixed lattice spacing we now take the results 
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L/a /3 



4 
4 
4 
6 
6 
6 



5.2000 
5.2000 
5.2000 
5.2000 
5.2000 
5.2000 



0.134700 
0.133780 

0.135600 
0.135200 



am 

-0.03745(41) 
-0.00086(35) 

-> 
-0.01322(26) 
+0.00289(24) 

-»• 



g 2 W 
3.730(11) 
3.797(11) 
3.798(11) 
4.810(32) 
4.984(33) 
4.954(33) 



Sit 



Li/a 



-0.686 5.11(3)(13) 



+0.470 



6 5.2638 
8 5.4689 
10 5.6190 



0.135673 
0.136575 
0.136700 



+0.00012(19) 
+0.00046(11) 
+0.00038(8) 



4.550(25) 
4.526(32) 
4.531(51) 



+0.066 
+0.042 
+0.037 



12 
16 



5.7580 
5.9631 



0.136623 
0.136422 



+0.00067(7) 
-0.00096(4) 



4.501(91) 
4.40(10) 



+0.017 
+0.084 



5.33(4)(11) 

5.89(4)(2) 
7.91(7)(1) 
9.87(14)(2) 

11.94(31)(1) 
16.40(50)(6) 



Table 13: Values of L\/a after correcting the simulated values L/a. For completeness we also 
include the two largest /3-values, where no new simulations enter. The second error on the final 
result is the systematic one. 

of g 2 (L) as computed at integer L/a and correct for the desired value L\/a at 
g 2 (Li) = 4.484, using 



In 



In 



+ 



(D.2) 



lg{ Ll ) ' 

The non-perturbative estimate of the /3-function from [5] is inserted, but we treat 
the difference to the value of Li/a, obtained from perturbation theory at the 
highest available order, as a systematic uncertainty. 10 We quote this as second 
error in our results for L\ja in Table 13, the first error is due to the statistical 
error of g 2 (L). Our final numbers to be used in Section 6.1 are those shown in 
the last column of Table 13. 



E Critical mass m cr 

Since Wilson fermions explicitly break chiral symmetry, the quark masses experi- 
ence an additive renormalization m cr . It is defined by the point where the PCAC 
quark mass eq. (3.3) vanishes for quarks degenerate to the sea quarks: myi = 0. 
Even though we have chosen our improvement procedure in terms of m rs and hence 
do not need m cr in our analysis, we give here its determination for completeness. 

10 The estimate of f)(x) in [5] applies to the continuum step scaling function, while we here 
apply it for not-so-large L/a. This represents a second reason to prefer L/a = 6 to L/a = 4 in 
our determination at f3 = 5.2. 
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Zr m 




5.2 


0.1360546(25) (30) 


1.438(20)(30) 


1.549(42) 


5.3 


0.1364572(11)(30) 


1.310(11)(20) 


1.323(23) 


5.5 


0.1367749(06) (04) 


1.228(06)(05) 


1.157(08) 



Table 14: Value of the critical hopping parameter and the renormalization constant Zr m and 
r m for our three values of the coupling constant. The first error in Zr m is statistical, the second 
systematic from a variation of the functional form and the range of the fit to eq. (E.l). For r m 
these errors are added in quadrature together with the one from Z. 

Up to terms of order a 2 the subtracted bare quark mass m q = m — m cr is 
proportional to the quark mass from the plateau of eq. (3.3) 

Zp 

m 12 (l + b T am 12 ) = Z m —r m m q = Zr m m Ci , (E.l) 

Za 

which we use to define the subtraction term m cr . The improvement coefficient b r is 
a linear combination of those improving the axial current, the pseudoscalar density 
and the bare subtracted quark mass while Z m renormalizes quark mass differences 
mq,i — Triqj and Z m r m renormalizes the trace of the bare subtracted mass matrix 
J2i m q_,i [30]. In one-loop perturbation theory, b r = | + 0.0500(?q [31,59], which we 
expect to be a sufficient approximation since our quark masses are small. Results 
are given in Table 14, where we quote k ct = (2am CT + 8) _1 instead of m cr and use 
the non-perturbative Z = 1 + 0.090514^ (1-0.3922 ^-0.2145 g$)/(l -0.6186^) 
of [46] to determine r m from Zr m . The central values are from fits of eq. (E.l) to 
the data with ami 2 < 0.01; the resulting x 2 /d.o.f. are all smaller than 1, which 
might be an indication that the errors on m are overestimated. The systematic 
error is estimated from including also a free cubic term on the right hand side of 
eq. (E.l) and a variation of the fit range up to amyi = 0.016. 

These non-perturbative values for Zr m and r m deviate significantly from the 
one-loop expectation Zr m = 1 + 0.090514^ [59], which is roughly Zr m ~ 1.1 
in our region of couplings and even more so for r m which is one up to 0((?q) 
corrections. 

F Masses of the pseudoscalar mesons, decay constants and 
quark masses 

Following the analysis described in Sect. 3.2, we determined starting time slices 
x™ m /a whose value for the ensembles we give in Table 15. We give one value each 
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for the mesons composed of two sea quarks, a sea and a strange quark and two 
strange quarks. We list the raw data of quark masses, pseudoscalar masses and 
decay constants in Tables 16-18. 
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24 


24 


26 
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31 


31 



Table 15: Values of x™ m /a used in the analysis, from top to bottom, for mesons made of two 
sea quarks, a sea and a strange quark and two strange quarks. 
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am r 



amps 



z^barc 
^PS 



A2 



0.13510 
0.13530 
0.13550 
0.13570 
0.13565 



0.13510 
0.13530 
0.13550 
0.13570 
0.13510 
0.13530 
0.13550 
0.13565 



0.02934 
0.02438 
0.01946 
0.01454 
0.02252 
0.02007 
0.01762 
0.01578 



0.3250 
0.2968 
0.2661 
0.2315 
0.2858 
0.2702 
0.2537 
0.2407 



0.1009(2) 
0.0977(2) 
0.0943(2) 
0.0908(2) 
0.0963(2) 
0.0947(2) 
0.0930(2) 
0.0917(2) 



A3 



0.13510 
0.13530 
0.13550 
0.13570 
0.13580 



0.13510 
0.13530 
0.13550 
0.13570 
0.13510 
0.13530 
0.13550 
0.13570 
0.13580 



0.02707 
0.02211 
0.01720 
0.01231 
0.01842 
0.01597 
0.01353 
0.01109 
0.00985 



0.3078 
0.2785 
0.2464 
0.2096 
0.2551 
0.2379 
0.2194 
0.1993 
0.1883 



0.0965(2) 
0.0930(2) 
0.0893(2) 
0.0854(2) 
0.0899(2) 
0.0882(2) 
0.0863(2) 
0.0843(2) 
0.0832(2) 



A4 



0.13510 
0.13530 
0.13550 
0.13570 
0.13590 



0.13510 
0.13530 
0.13550 
0.13570 
0.13510 
0.13530 
0.13550 
0.13570 
0.13590 



0.02557 
0.02061 
0.01571 
0.01087 
0.01574 
0.01330 
0.01087 
0.00845 
0.00601 



0.2958 
0.2659 
0.2328 
0.1948 
0.2334 
0.2150 
0.1949 
0.1726 
0.1466 



0.0933(2) 
0.0896(2) 
0.0857(2) 
0.0814(2) 
0.0851(2) 
0.0833(2) 
0.0813(2) 
0.0791(2) 
0.0766(3) 



A5 



0.13510 
0.13530 
0.13550 
0.13570 
0.13594 



0.13510 
0.13530 
0.13550 
0.13570 
0.13510 
0.13530 
0.13550 
0.13570 
0.13594 



0.02488 
0.01993 
0.01504 
0.01021 
0.01460 
0.01217 
0.00975 
0.00734 
0.00444 



0.2907 
0.2605 
0.2269 
0.1882 
0.2244 
0.2053 
0.1844 
0.1608 
0.1263 



0.0918(2) 
0.0881(2) 
0.0841(2) 
0.0797(2) 
0.0828(2) 
0.0810(2) 
0.0789(2) 
0.0766(3) 
0.0735(3) 



Table 16: Partially quenched average quark masses, pseudoscalar meson masses and decay 
constants in lattice units for (3 = 5.2. 
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am rs 


amps 




7-iKa rp 

^PS 




0.13540 





13540 


0.03121(7) 


0.3018 


;i2) 


0.0884(5) 




0.13560 





13560 


0.02589(8) 


0.2746 


13) 


0.0852(5) 




0.13580 





13580 


0.02061(8) 


0.2449 


;i3) 


0.0819(5) 


E4 


0.13600 





13600 


0.01535(8) 


0.2118 


,14) 


n nTo A / f \ 

0.0784(5) 




0.13610 





13540 


0.02191(8) 


0.2529 


,14) 


n no r\ r~ 1 f \ 

0.0825(5) 









13560 


r\ Ai n O O / O \ 

0.01928(8) 


0.2372 


14) 


n no nn { r* \ 

0.0809(5) 









13580 


0.01666(8) 


0.2206 


,14) 


0.0793(5) 









13600 


0.01404(9) 


0.2028 


,15) 


0.0775(5) 









13610 


0.01272(9) 


0.1934 


,15) 


r\ r\'-7 r"-7 ( n \ 

0.0767(6) 




0.13540 





13540 


0.02968(3) 


0.2896 


;3) 


0.08509(17) 




0.13560 





13560 


0.02436(3) 


0.2619 


;3) 


0.08173(17) 




0.13580 





13580 


0.01907(3) 


0.2318 


;3) 


0.07809(17) 


E5 


0.13600 





13600 


0.01383(3) 


0.1979 


' A \ 

,4) 


0.07413(17) 




0.13625 





13540 


0.01840(3) 


0.2281 


'4) 


0.07672(17) 









13560 


r\ r\-\ r ^7n / o\ 

0.01578(3) 


0.2114 


'4) 


0.07514(17) 









13580 


0.01317(3) 


0.1934 


' A \ 

,4) 


0.07337(17) 









13600 


0.01056(3) 


0.1737 


'4) 


0.07139(17) 









13625 


0.00727(3) 


n 1 A V A 

0.1454 


' A \ 

,4) 


0.06852(20) 




0.13540 





13540 


0.02874(2) 


0.2824 


;s) 


0.08294(18) 




0.13560 





13560 


0.02341(2) 


0.2542 


'3) 


0.07947(18) 




0.13580 





13580 


0.01812(3) 


0.2234 


[3) 


0.07570(19) 


F6 


0.13600 





13600 


f~\ A1 00 0/ 0\ 

0.01288(3) 


n i o o r* 

0.1885 


,3) 


0.0716(2) 




0.13635 





13540 


0.01616(3) 


0.2113 


'4) 


0.0729(2) 









13560 


0.01354(3) 


0.1935 


^4) 


0.0713(2) 









13580 


0.01092(3) 


0.1741 


,4) 


0.0695(2) 









13600 


r\ AAAfif\/n\ 

0.00832(3) 


0.1524 


'4) 


0.0674(2) 









13635 


0.00374(3) 


0.1036 


'5) 


0.0629(3) 




0.13540 





13540 


0.028455(17) 


0.2801 


;2) 


0.08208(14) 




0.13560 





13560 


0.023121(19) 


0.2520 


;3) 


0.07857(14) 




0.13580 





13580 


0.017832(19) 


0.2211 


;3) 


0.07473(15) 


F7 


0.13600 





13600 


0.012594(20) 


0.1860 


;3) 


0.07045(16) 




0.13638 





13540 


0.01549(2) 


0.2068 


;3) 


0.07122(18) 









13560 


0.01288(2) 


0.1886 


;3) 


0.06957(18) 









13580 


0.01027(2) 


0.1687 


;3) 


0.06771(18) 









13600 


0.00767(2) 


0.1463 


;3) 


0.06557(18) 









13638 


0.00272(2) 


0.0886 


[A) 


0.0603(2) 



Table 17: Partially quenched average quark masses, pseudoscalar meson masses and decay 
constants in lattice units for (3 — 5.3. 
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am rs 


amps 




7 — inavr 1 




0.13600 





13600 


0.02308(3) 


0.2172 




0.0631(4) 




0.13615 





13615 


0.01885(3) 


0.1958 


y) 


0.0605(4) 




0.13630 





13630 


0.01462(3) 


0.1723 


{&) 


0.0578(4) 


N4 


0.13645 
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0.01040(3) 


0.1456 


8) 


r\ r\ r A ( A \ 

0.0548(4) 




0.13650 





13600 


0.01602(3) 


0.1805 


'%) 


0.0584(4) 









13615 


0.01391(3) 


0.1682 


J) 


0.0572(4) 









13630 


0.01181(3) 


0.1550 


'8) 


0.0558(4) 









13645 


0.00970(3) 


0.1407 


'9) 


0.0543(4) 









13650 


0.00899(3) 


0.1357 


8) 


0.0538(4) 




0.13600 





13600 


0.02265(3) 


0.2144 


7) 


0.0622(4) 




0.13615 





13615 


0.01841(3) 


0.1927 


y) 


0.0595(4) 




0.13630 





13630 


0.01419(3) 


0.1689 


y) 


0.0566(4) 


N5 


0.13645 





13645 


0.00997(4) 


0.1417 


'8) 


0.0534(4) 




0.13660 





13600 


0.01417(4) 


0.1689 


7) 


0.0558(4) 









13615 


0.01207(4) 


0.1559 


y) 


0.0546(3) 









13630 


0.00996(4) 


0.1418 


y) 


0.0532(3) 









13645 


0.00786(4) 


0.1262 


'8) 


0.0516(3) 









13660 


0.00575(4) 


0.1085 


'8) 


r\ r\ A r\c\ / A \ 

0.0499(4) 




0.13600 





13600 


0.022299(17) 


0.2099 




0.0609(3) 




0.13615 





13615 


0.018062(17) 


0.1881 




0.0581(3) 




0.13630 





13630 


0.013835(18) 


0.1642 




0.0552(2) 


N6 


0.13645 





13645 


0.009619(19) 


0.1369 


[6) 


0.0519(3) 




0.13667 





13600 


0.01283(2) 


0.1586 


'6) 


0.0534(3) 









13615 


0.01073(2) 


0.1450 


'5) 


0.0522(3) 









13630 


0.00863(2) 


0.1301 


^5) 


0.0507(3) 









13645 


0.00653(2) 


0.1135 


^5) 


0.0491(3) 









13667 


0.00343(3) 


0.0834 


y) 


0.0461(3) 




0.13600 





13600 


0.022117(10) 


0.2085 


;3) 


0.06040(14) 




0.13620 





13620 


0.016468(10) 


0.1789 


'3) 


0.05666(14) 




0.13640 





13640 


0.010837(12) 


0.1447 


v 3 ) 


0.05241(14) 


07 


0.13660 





13660 


0.005226(14) 


0.1012 




0.04741(15) 




0.13671 





13600 


0.012084(13) 


0.1535 


;3) 


0.05225(18) 









13620 


0.009283(13) 


0.1345 


;3) 


0.05050(16) 









13640 


0.006483(13) 


0.1126 


;3) 


0.04841(16) 









13660 


0.003681(15) 


0.0855 


'3) 


0.04581(16) 









13671 


0.002130(15) 


0.0658 


[A] 


0.0440(2) 



Table 18: Partially quenched average quark masses, pseudoscalar meson masses and decay 
constants in lattice units for (3 — 5.5. 
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